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ABSTRACT 

We present previously unpublished July 2005 iJ-band coronagraphic data of the young, planet- 
hosting star HR 8799 from the newly-released Keck/NIRC2 archive. Despite poor observing condi- 
tions, we detect three of the planets (HR 8799 bed), two of them (HR 8799 be) without advanced 
image processing. Comparing these data with previously published 1998-2011 astrometry and that 
from re-reduced October 2010 Keck data constrains the orbits of the planets. Analyzing the planets' 
astrometry separately, HR 8799 d's orbit is likely inclined at least 25° from face-on and the others may 
be on in inclined orbits. For semimajor axis ratios consistent with a 4:2:1 mean-motion resonance, 
our analysis yields precise values for HR 8799 bed's orbital parameters and strictly constrains the 
planets' eccentricities to be less than 0.18-0.3. However, we find no acceptable orbital solutions with 
this resonance that place the planets in face-on orbits; HR 8799 d shows the largest deviation from 
such orbits. Moreover, few orbits make HR 8799 d coplanar with b and c, whereas dynamical stability 
analyses used to constrain the planets' masses typically assume coplanar and/or face-on orbits. This 
paper illustrates the significant science gain enabled with the release of the NIRC2 archive. 
Subject headings: planetary systems, stars: early-type, stars: individual: HR 8799 



1. INTRODUCTION 

The nearby, yo ung A- type star HR 8 799 (d = 39.4 
pc, ps 30 Myr; iZuckerman et all 1201 If) harbors the 
first independently confirmed, directly imaged exoplan- 
etary system and t he only imaged multi-planet system 
(jMarois et al.1 120081 ) . After the discovery of HR 8799 
bed (rnmi 24, 38, and 68 AU) reported in Novem- 
ber 2008 (|Marois et al.ll2008f ). other studies identified at 
least one of these planets in ar c hival data taken prior 
to 2008 (iLafreniere et al I I2009t iFukagawa etail 120091: 
iMetchev et al.ll2009HSoummer et al.ll201lD . 

HR 8799 planet astrometry derived from both pre 
and post-discovery images can help constrain the sys- 
tem's dynamical stability and, in turn, the planets' phys- 
ical properties. At least two of the HR 8799 plan- 
ets are likely locked in a mean motion resonance, oth- 
erwise th e system would quickly become dynamically 
unstable (jFabrvckv and Murrav-Clavl l2010f ). The re- 
cently discovered fourth companion at ~ 15 AU, HR 
8799 e, generally makes dynamical stab ility less likely 
(jMarois et al.ll2010at ICurrie et al.ll2011al) . favoring lower 
masses of Mb jC dc < 7, 10 Mj, an important constraint 
given the uncertainties in derivin g masses from planet 
cooling and atmosphere mode ls (jSpiegel and Burrowsl 
[2011 iMadhusudhan et al1l27)TTh . 

Studies focused on fitting the planets' orbits and/or 
testing dynamical stability typically assume that the 
planets are a) in resonance (4:2:1 for HR 8799 bed 
or 2:1 for HR 8799 cd), b) in cir cular, face-on orbits , 
c) and/or in coplanar orbits (e.g. iMarois et al.l l2010at 
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ICurrie et al.l l2011al see also Fabrycky a nd M urray- 
Clay 2010). However, iSoummer et all (|2011h show 
that circular, face-on, and coplanar orbits are incon- 
sistent with 1998 HST astrometry, identifying a best- 
fit orbit for HR 8799 d of i = 28° and e = 0.115. 
Generally, more eccentric orbits destabilize the system. 
The system stability depends on the (mutual) inclina- 
tions of the planets (iFabrvckv and Murrav-Clavl 120101 : 
iSudol and Haghighipourll2012D . Thus, the HR 8799 plan- 
ets' true mass limits derived from dynamical stability 
arguments may slightly differ from those previously re- 
ported. 

Well-sampled HR 8799 d astrometry could help clarify 
whether HR 8799 d's orbit must be inclined, eccentric, 
and/or coplanar with the other planets. However, until 
now there is a ~ 9-year gap between the 1998 HST detec- 
tion and the next one (2007; Metchev et al. 2009). New 
astrometry for HR 8799 bee in between 1998 and 2007 
could also help constrain those planets' orbits. By bet- 
ter determining the HR 8799 planets' orbital properties, 
we can more conclusively investigate system dynamical 
stability and thus better clarify the range of allowable 
planet masses. 

In this Letter, we report the detection of HR 8799 bed 
from unpublished, now-public Keck/NIRC2 data taken 
in 2005 supplemented w ith a re-reduc t ion of published 
October 2010 data from IMarois et al.l (|2010aD . We use 
these data to better constrain the orbital properties of 
HR 8799 bed. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. July 2005 Data 

We downloaded HR 8799 data taken on July 15, 
2005 from the newly- available Keck/NIRC2 data archive 
(Program ID H53BN2, P.I. Michael Liu). The data 
were taken in H band with the narrow camera (9.952 
mas /pixel: lYelda et al.ll2010h with the 0'.'6 diameter coro- 
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nagraphic spot and the "incircle" pupil plane mask. HR 
8799 was observed in 10 s exposures in "vertical angle" 
or angular differential imaging mode (|Marois et al.ll2006ft 
through transit (Hour Angle = [—0.20,0.55]) with a total 
field rotation of 147.1°. During these observations, the 
seeing conditions fluctuated and the observers periodi- 
cally recentered the star behind the mask, changing the 
intensity profile of the stellar halo (and thus the quasi- 
static speckle pattern) . During a few frames near transit 
the star did not properly center at all behind the coron- 
agraph. We identify ~ 11 minutes of science-grade data. 
Basic image processing followed s tandard steps pre vi- 
ously used to process NIRC2 data (jCurrie et alj |2012). 

For a first-order reduction, we perform a simple "clas- 
sical" ADI-based PSF subtraction (e.g. iMarois et al.l 

2006) . Figure Q] shows this reduction, clearly reveal- 
ing HR 8799 b (SNR - 12) and identifying HR 8799 
c, albeit at low SNR (~ 4). With the LOCI approach 
(|Lafreniere et al.ll2007D as implemented and modified in 
previous work (|Currie et al.|[20Toll201 lalfbh ■ we easily de- 
tect HR 8799 b and c and obtain a marginal detection of 
HR 8799 d at r ~ 0'.'6 several degrees a way (clockwise) 
from t he July 2007 position reported by IMetchev et al.l 
(2009) (not shown). 

To improve the signal-to-noise of our HR 8799 d detec- 
tion, we incorporated several upgrades to enhance con- 
trast which are being implemented in a new "adaptive" 
LOCI (A-LOCI) pipeline (T. Currie, 2012 in prep.; see 
also similar steps in Marois et al. 2010b). We subtract off 
the seeing halo in each image to measure the static/quasi- 
static speckle pattern, determining the cross-correlation 
function for the speckle patterns in annular sections for 
all possible image pairs. While the speckle pattern is 
generally better correlated between frames taken close 
together in time, this is not always the case, especially 
when comparing frames before and after small telescope 
nods. Therefore, for each annular section of a science 
image we want to subtract, we filter reference image sec- 
tions by their degree of correlation to remove those below 
a certain, predetermined threshold (r corr ). 

HR 8799 b is detectable in most individual processed 
frames (SNR/frame ~ 4-7), so we measure its position 
to identify and correct for any astrometric biases caused 
by a PA "jump" for frames obtained near transit due to 
imperfect mechanical alignment of the telescope's y axis. 
The position angles of HR 8799 b in frames more than 
0.25 hours from transit are consistent, but the PA offset 
follows a bell-shaped curve with a maximum offset of ks 
0.6 degrees centered on transit. We model and correct 
for this offset using a fifth-order polynomial. We also 
reran our pipeline with different rotation axis offsets due 
to image registration errors, setting an upper limit to 
this of 0.5 pixels in each coordinate. 

Furthermore, unlike the original LOCI algorithm, we 
set the azimuthal length of the subtraction zone to be 
smaller (not identical) to the azimuthal length for the op- 
timization zone, equal to dr (the length along the radial 
direction). We then center the optimization zone on the 
subtraction zone. Finally, we iterativel y determine the 
algor ithm parameters - 5, Na, <?, dr ( Lafreniere et alJ 

2007) , and r corr - that maximize the signal-to- noise of 
point sources and applied these settings to extract our 



final imagefl 

Figure [T] (right panel) shows our final image displaying 
higher signal-to-noise detections of HR 8799 b (SNR ~ 
38) and c (SNR ~ 18) and new detection of d (SNR - 5), 
using algorithm parameters of 5 > 0.74, Na — 245, g = 
0.95, dr = 6, and r corr > 0.315, though similar settings 
yield nearly identical results (i.e. 5 > 0.73-0.8, r corr > 
0.31-0.36). We achieve contrast gains of up to ~ 80% 
(for HR 8799 d) over our best LOCI reduction. 

For flux calibration, we perform aperture photome- 
try on HR 8799 bed and on the nearby star GJ 616.2 
observed just prior to HR 8799. We use fake point 
sources to further correct for LOCI-based photomet ric 
biases (i.e. lLafreniere et~aLl [20071: iCurrie et al.ll2011allbl ) 
and find m(H) = 18.05 ± 0.09 mag for HR 8799 b, m{H) 
= 17.06 ± 0.13 for HR 8799 c, and m(H) = 16.71 ± 0.24 
for HR 8799 d. The magnitude differences between HR 
8799 c and d appear slightly discrepant co mpared to bet- 
ter ca librated NIRC2 measurements from IMarois et al.l 
(2008), although the individual measurements are con- 
sistent to within ~ 1 a. Photometry derived for HR 8799 
be using classical PSF subtraction agrees with that de- 
rived from our A-LOCI based reduction within errors. 
Bright residual speckles at r = / /3-0 / /4 prevent detect- 
ing HR 8799 e. 

2.2. October 2010 Data 

To supplement the 2005 HR 8799 astrometry, we down- 
loaded and reduced October 2010 L'-band NIRC2 data 
from the Keck archive (P. I. B. Macinto s h). Th ese data 
are the latest reported by IMarois et al.l ()2010al ) who fo- 
cus on HR 8799 e astrometry: we use these data in- 
stead to extract astrometry for HR 8799 bed. Individual 
exposures consist of 50 s frames totaling ps 80 minutes 
taken through transit without a coronagraph with peri- 
odic telescope nods for sky subtraction. Observing con- 
ditions appeared variable at a level comparable to the 
July 2005 data and worse than other recen t Keck data 
(i.e. ICurrie etai]|20T3 iRodigas et al.l[2012h . but we de- 
tect all four planets (SNR ~ 6-20) with our pipeline. 

3. ASTROMETRIC ANALYSIS AND ORBIT FITTING 
3.1. Method 

We use our detections to better constrain the orbits 
of HR 8799 bed, first by fitting the orbits of the planets 
separately, then identifying the subset of orbits consis- 
tent with systematically more stable mean motion res- 
onance configuration. We calibrate our astrometry by 
assessing and correcting for biases introduced by LOCI- 
based processing in the same manner as our photomet- 
ric calibration, using fake point sources. Using different 
telescopes and slightly different image processing tech- 
niques lead to systematic biases in planet astrometry. 
To minimize these biases, we restrict ourselves to as- 
trometry from HST /1998 (iSoummer et all 120111), Sub- 
aru/20 02 and 2009 (|Fukagawa et alJ l2009t ICurrie et al 



l2011ah . Keck/2004, 2005, 2007-2009 (|Marois et al 



6 We also considered a "reference PSF library" derived from 
other 2004-2005 NIRC2 data obtained with the same setup (coro- 
nagraph size, filter, etc.) to further attenuate speckles. However, 
this library degrades the SNR of HR 8799 b and c by ^ 30-60% 
and renders HR 8799 d undetectable, because HR 8799's speckle 
patterns are poorly correlated with library's. 
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[20081 IMetchev et all [200l IGalicher et alJ [20Til this 
work), VLT/ 200 9 (ICurrie et al.l l2011aTlmd LBT/Pisces 
(|Esposito et al.l I2012I ). For the Pisces data, we in- 
clude the substantial north position angle uncertainty. 
We furth ermore modify the Keck/NIRC2 astrome try re - 
ported in lMarois etail (|2008[ ) and lMetchev et alJ (2009) 
to reflect the upd ated NIRC2 astrometric calibration 
(jYelda et alJ 12010) accuracy, rescaling the position by 
a factor of 9.952/9.963, and putting in a PA nort h,ncvj- 
-P^north.oid = 0.13° clockwise rotation. 

To separately determine the range of allowable HR 
8799 bed orbits, we follow our previously-adopted 
Monte Carlo-based approach (Thalma nn et al.1 120091 : 
ICurrie et al. l 1201 lbf) . comparing the HR 8799 planet as- 
trometry to predictions from randomly-selected orbits. 
In a first set of "conservative" simulations, we con- 
sider the orbits separately. Dynamical stability analy- 
sis suggests that HR 8799 bed a re likely in resonance 
(|Fabrvckv and Murray- Clavl |2010| ) : the 4:2:1 resonance 
is particularly adept at stabilizing the system. To focus 
on dynamically stable orbits we then select the subset of 
best-fitting orbits that preclude the planets from cross- 
ing orbits and are consistent with a 4:2:1 resonance. Here 
we define "resonance" broadly, including orbits with ra- 
tios of periods between 1.9 and 2.1 for consecutive pairs 
of planets since, at least in s ome circumstances, e xact 
period ratios may be rare (e.g. iFabrvckv et alll^Oll) Q. 

For all our simulations, the minimum \ 2 for the b, 
c, and d planets are 20.8, 13.8, and 14.8, for reduced 
\l valu es of 1.04, 0.86, and 1.24. Following ICurrie et al.1 
(l2oTiB . we choose a cutoff of \ 2 < Xmin + 1 to represent 
the family of best-fitting orbits. Formally, this cutoff 
admits only an average additional deviation per each x 
or y measurement of « l/2iV obs (i.e. w 1/20 for HR 8799 
b; » 1/14 for HR 8799 d) beyond the best-fit models 
which themselves imply typical deviations of ~ 1-a per 
each x or y measurement (since the minimum reduced 
xt values are « 1). However, we obtain nearly identical 
results for more relaxed cutoffs (see below) . From the set 
of models passing our x 2 cutoff, including the subset in 
resonance, we determine the weighted median value and 
the weighted 68% confidence interval about the median 
for each model parameter from amongst the set of best- 
fitting orbits. 

Furthermore, we report a "most likely orbit" (MLO) 
simply defined as follows. First, over the best-fitting 
family of orbits, we calculate histograms for the follow- 
ing parameters: logarithm of the semi-major axis (h\ oga ), 
eccentricity (h e ), inclination (hi), longitude of ascending 
node (/in), and argument of periastron (h^). For each 
orbit n in the best-fitting ensemble, we then define the 
measure of likelihood C(n) as 

C(n) = hi oga (loga n ) ■ h e (e n ) ■ hi(i n ) ■ ha(fl n ) -h^ojn) ■ W, 

i.e., the product of all histograms values in the bins in 
which the orbit n lies, representing the individual like- 
lihood of each measured orbital parameter, as well as 
the statistical weight, W, representing the likelihood of 

7 Dynamical simulations identify stable solutions for at least 
some planet masses where only HR 8799 c and d are in resonance, 
so our conclusions from this set of astrometric analyses may be less 
applicable for HR 8799 b. 



the observed planet position within the orbit (i.e., the 
anomaly j3- The most likely orbit is then the one orbit 
that maximizes the measure of likelihood, £(jimlo) = 
max n C(n). Because of the highly skewed distribution of 
some parameters from the best-fit orbits, in particular 
log a, the MLO parameters can differ significantly from 
the weighted median parameters. 

3.2. Results 

Table 1 summarizes our results, and Figure [5] displays 
the orbits in a/i/e space (left), the histogram distribu- 
tion of i (middle), and the histogram distribution of the 
longitude of the ascending node, fi (right). The top pan- 
els display properties for HR 8799 b, the middle for HR 
8799 c, and the bottom for HR 8799 d. Assuming a x 2 
cutoff of Xmin + 1) the observations well constrain the 
HR 8799 d orbit - a ~ 24-32 AU, i ~ 32-42°, e ~ 0.03- 
0.23, and fl ~ 43-63°. They limit the HR 8799 c's most 
plausible orbital parameters to a ~ 36-42 AU, i ~ 13- 
26°, and e ~ 0.03-0.13. As expected, the parameters for 
HR 8799 b are the most poorly constrained, showing the 
widest dispersion and the largest differences between the 
median parameter value and that from the MLO. 

Our analysis clearly disfavors face-on orbits for all 
three planets, especially for HR 8799 c and d. Addition- 
ally, the inclination distributions for HR 8799 be appear 
systematically skewed towards values lower than those 
for d. Formally, though, the set of acceptably-fitting or- 
bits HR 8799 bed include some that make the planets 
coplanar. 

Assuming a 4:2:1 mean-motion resonance configura- 
tion, we place far stronger limits on nearly all of the HR 
8799 bed orbital properties (Figure [3]). This assumption 
explicitly rules out e > 0.18 for HR 8799 be and e > 0.3 
for HR 8799 d. Likewise, we identify a very narrow range 
of planet semimajor axes: a = 67.5-70.8 AU, 42.1-44.4 
AU, and 26.4-28.1 AU for HR 8799 b, c, and d. The 68% 
confidence interval in Q for HR 8799 d further narrows 
to 46°-62°. 

For this configuration, HR 8799 d (c) must be in an 
orbit viewed more than 25° (15°) from face-on while HR 
8799 b is likely inclined by at least 5°. Furthermore, 
the inclination distributions between HR 8799 d and HR 
8799 be are even more dissimilar, implying that HR 8799 
d is most likely inclined relative to c by at least rs 7° 
and b by more than ss 21°. While our analyses cannot 
conclusively rule out coplanar orbits in a 4:2:1 mean- 
motion resonance, they suggest that few such orbits are 
compatible with 12 years of HR 8799 planet astrometry. 

To confirm that we are fully sampling the subset of 
orbital parameters covering the x 2 minima, we run our 
simulations with a more relaxed x 2 cutoff of % 2 < Xmin + 
[5, 5, 3.5] for HR 8799 b, c, and d, which formally admit 
an additional average deviation from the data of « 0.3- 
a for each measurement. With this cutoff, we obtain 
nearly identical results (second set of rows in Table 1). 
Considering the planets' orbits separately, we find a ~ 
23-31 AU, i ~ 27-41°, e - 0.01-0.31, and Q ~ 41-53° for 
HR 8799d. We find similar ranges in orbital parameters 
for HR 8799 c and (for the resonance case) HR 8799 

8 Here, the statistical weight W is denned as the mean orbital 
velocity for the corresponding orbit divided by the orbital velocity 
at the observed epoch, W := (^)orbit/' u (iobs)- 
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b. Likewise, the planets' range of inclinations exclude 
face-on orbits. HR 8799 d's inclination distribution is 
skewed to systematically higher values, expected if it is 
non-coplanar with b and c, although here there are more 
orbit combinations that could make the planets coplanar. 

4. DISCUSSION 

From analyzing HR 8799 bed astrometry from our new 
"pre-discovery" image and other data, we provide new 
constraints on the planets' orbital properties. Treating 
the three planets separately, we narrowly constrain three 
major orbital parameters (a/i/e) for HR 8799 c and d. 
None of the planets are likely to be orbiting face-on and 
the inclinations for acceptably-fitting orbits are system- 
atically higher for HR 8799 d than for HR 8799 b and 
c. 

If HR 8799 bed have semimajor axes consistent with 
a 4:2:1 resonance, our analysis strongly constrains the 
major orbital properties for all three planets. The three 
planets (especially c and d) then even more obviously 
have inclined orbits. Most acceptable solutions for HR 
8799 d place the planet on an orbit inclined by more than 
7° (21°) relative to HR 8799 b(c)'s orbit: few orbital so- 
lutions consistent with the astrometry also place them 
on coplanar orbits. Adopting a less restrictive defini- 
tion for "acceptably-fitting" orbits does not undo any of 
these trends, although there are more orbit combinations 
making the planets coplanar. Adopting the median pa- 
rameter value or MLO instead of the more conservative 
68% confidence interval likewise does not change these 
results. 

These results provide valuable input for constraining 
the mass of the HR 8799 planetary system. Longer-term 
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astro metric monitoring of HR 8799 (i.e. lKonopackv et al.l 
1201 If ) will better clarify the planets' orbital properties. 
Limits on the planets' dynamical masses will provide cru- 
cial input for planet cooling models and even more firmly 
establish HR 8799 as a benchmark system to understand 
the properties of young, self-luminous planets. 

Finally, this work and other recent studi es of 

HR 8799 (iSoummer et all [20Tll: lLafreniere et~aTl [2Q0l 
iFukagawa et al.l 12009ft clearly demonstrate the value of 
publicly archiving data on advanced telescopes. In our 
case, detecting at least two HR 8799 planets (HR 8799 
be) was rather straightforward and did not require ad- 
vanced image processing techniques developed well after 
the data were taken. As data for Keck and many other 
8-10 m class telescopes are now archived, they provide an 
indispensible resource with which to confirm and charac- 
terize directly imaged planets like HR 8799's and other 
substcllar companions. 
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TABLE 1 

HR 8799 Planet Photometric/ Astrometric Properties 



Parameter 


HR 8799 b 


HR 8799 c 


HR 8799 d 


Measured 

Photometry 
m(H) on 2005-07-15 
Astrometry 
2005-07-15 ([E,N])" 
2010-10-30 ([E,N])" 


18.05± 0.09 

1.496, 0.856 (± 0.005) 
1.546, 0.748 (± 0.005) 


17.06 ± 0.13 

-0.713, 0.630 (± 0.005) 
-0.598, 0.737 (± 0.005) 


16.71 ± 0.24 

-0.087, -0.578 (± 0.010) 
-0.283, -0.567 (± 0.005) 



Derived 

MLO, med.,[68% C.I.] 

(Xfijn < X^ in + 1) 

a (AU), full 

// 4:2:1 resonance 

i (°), full 

// 4:2:1 resonance 
e, full 

// 4:2:1 resonance 

n (°), mi 

II 4:2:1 resonance 



(xf im < x 2 



*lim ^ A min -r[5,5,3.5]) 

a (AU), full 
// 4:2:1 resonance 
i (°), full 
// 4:2:1 resonance 
e, full 

// 4:2:1 resonance 

n (°), mi 

II 4:2:1 resonance 



71.0, 109.9 [69.7,164.9] 
68.1, 68.8 [67.5,70.8] 
14.1, 34.9 [12.3,43.0] 

8.5, 9.5 [4.9,14.8] 
0.02, 0.27 [0.02,0.49] 
0.01, 0.02 [0.0,0.03] 
149.7, 141.0 [40.9,161.4] 
163.2, 87.1 [22.9,158.9] 



70.7, 80.8 [68.2,117.2] 
68.0, 68.5 [66.4,71.0] 
15.2, 24.5 [10.4,36.6] 
10.8, 11.5 [6.1,17.1] 
0.01, 0.13 [0.02,0.34] 
0.01, 0.03 [0.01,0.06] 
142.7, 119.5 [32.4,160.5] 
142.1, 78.2 [21.9,154.3] 



37.2, 38.0 
42.5, 43.2 

21.7, 19.9 

25.8, 27.6 
0.01, 0.08 
0.12, 0.14 
128.8 
147.4 



35.5,42.0 
42.1,44.4 
12.5,25.8' 
25.2,28.8 
0.03,0.13 
0.11,0.17 
122/2 [62.0,152.4] 
131.9, [104.8,158.5] 



38.2, 39.8 
42.7, 42.9 

19.7, 23.2 

28.8, 27.4 
0.01, 0.07 
0.02, 0.09 
124.7, 86.8 [37. 
61.1, 70.5 [42.2 



36.6 
41.3 
14.0 
24.6 
0.02 
0.02 



,46.2 
44.6 
,31.3 
.30.2 
,0.15 
,0.14 
1,142.5] 
,140.5] 



26.2, 27.3 

27.3, 27.3 
37.1, 37.9 
38.1, 37.9 
0.04, 0.09 
0.03, 0.04 

56.8, 53.5 

56.9, 54.3 



26.2, 27.0 

27.3, 27.1 
38.1, 36.9 
37.1, 37.7 
0.01, 0.15 



24.4,31.5 
26.4,28.1 
31.6,41.6' 
36.0,39.1 
0.03,0.23 
0.01,0.08 
43.3,63.0 
46.0,60.2 



23.0,31.0 
26.0,28.2 
27.1,41.2' 
34.6,39.9' 
0.04,0.32 



0.0, 0.05 [0.01,0.11] 
52.5, 58.4 [41.0,98.5] 
60.9, 55.7 [44.1,67.2] 



Note. — Measured parameters — Our photometric uncertainties consider both the signal-to-noise of our 
detections and the absolute flux calibration uncertainties; astrometric uncertainties consider the SNR, astro- 
metric calibration uncertainty (e.g. 0.5 pixels in x and y, see Sect. 2), etc. Derived parameters — The three 
column entries are MLO (the "most likely orbit" (see §3), med. (the median parameter value), and 
C.I.] (the 68% confidence interval). 





Fig. 1. — HR 8799 images processed with "classical" PSF subtraction (left) and A-LOCI (right) showing the detections of HR 8799 b, 
c, and d (circled). 
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Fig. 2. — Astrometric analysis results for HR 8799 b (top), HR 8799 c (middle), and HR 8799 d (bottom) considering the full range of 
orbits satisfying the criterion xi < xt rain ^ ne P ane l s show the orbits in a/e/i space (left) and histogram distributions of the orbital 
inclination i (middle panels) and longitude of ascending node (right panels). In the left panels, the 'x' identifies parameters e and a p from 
the best-fit orbital solution; the 'plus' sign denotes the weighted median value for the same parameters. For the middle and right panels, 
the vertical blue dashed line identifies the MLO, the vertical red dashed line identifies the median parameter value, and the vertical red 
dotted lines bracket the 68% confidence interval. 
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Fig. 3. — Same as Figure [2] except only for non-crossi ng orbits consistent with a 4:2:1 m ean-motion resonance between HR 8799 b, c, and 
d, a configuration which promotes orbital stability (c.g Fabrycky and Murrav-Clav 2010). Note the lack of high-eccentricity orbits and the 
narrower range in acceptable orbital parameters, especially for HR 8799 c and d. 



